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In the linear response approximation, coherently driven microcavity polaritons in the pump-only 
configuration are expected to satisfy the Landau criterion for superfluidity at either strong enough 
pump powers or small flow velocities. Here, we solve non-perturbatively the time dependent Gross- 
Pitaevskii equation describing the resonantly-driven polariton system. We show that, even in the 
limit of asymptotically large densities, where in linear response approximation the system satisfies 
the Landau criterion, the fluid always experiences a residual drag force when flowing through the 
defect. We illustrate the result in terms of the polariton lifetime being finite, finding that the 
equilibrium limit of zero drag can only be recovered in the case of perfect microcavities. In general, 
both the drag force exerted by the defect on the fluid, as well as the height of Cherenkov radiation, 
and the percentage of particles scattered by the defect, show a smooth crossover rather than a sharp 
threshold-like behaviour typical of superfluids which obey the Landau criterion. 



PACS numbers: 03.75.Kk, 71.36+c, 41.60.Bq 

I. INTRODUCTION 

In the past two decades, microcavity polaritons^ have 
attracted considerable interest because of the possibil- 
ity of strongly coupling light and matter, leading to an 
easy manipulation and detection of the non-linear prop- 
erties of matter via light — see, e.g., RefPand references 
therein. Recently, in light of their out-of-cquilibrium na- 
ture, there has been a growing interest in studying polari- 
ton superfluid propertied. In fact, polaritons have very 
short lifetimes (around ~ 30 ps even in the best avail- 
able samples^) , therefore any polariton fluid is the result 
of a steady state balance between pumping and decay. 
The effects of pump and decay in polariton systems are 
the subject of several recent theoretical and experimental 
work. 

For condensates in local thermal equilibrium, such 
as superfluid 4 He and the ultracold atom Bose-Einstein 
condensates, the concept of superfluidity is strongly 
linked to several paradigmatical properties, such as the 
Landau criterion, quantised circulation of velocity, and 
metastable persistent flow. In particular, the Landau 
criterion connects the frictionless motion of a defect at 
velocities smaller than a critical one, with the shape of 
the spectrum of elementary excitations^. For weakly in- 
teracting Bose gases and a microscopic weak defect, such 
that a perturbative linear response theory can be ap- 
plied, the phonon-like dispersion of the Bogoliubov spec- 
trum demands that the critical velocity for dissipationless 
flow coincides with the velocity of souncP. However, for 
macroscopic defects, the critical velocity for the onset of 
drag is smaller than the speed of sound and most likely 
related to vortex nucleatiorP. 

In contrast, the relevance of the Landau criterion 
for out-of-cquilibrium condensates is questionable, not 
least because the spectrum of excitations is now complex 
rather than real. For polariton fluids, one has to singu- 
larly assess the system properties in the three different 



pumping schemes available: i) non-resonant pumping; 
ii) parametric drive in the optical-parametric-oscillator 
(OPO) regime; hi) coherent drive in the pump-only con- 
figuration. 

As far as the spectrum of quasi-particle excitations 
and the Landau criterion are concerned, the cases of 
non-resonantly pumped polaritons^ and parametrically 
driven polaritons in the OPO regime^ are similar: in both 
cases, there is a U(l) phase symmetry which is spon- 
taneously broken above a pump power threshold. This 
leads to the appearance of a gapless (Goldstone) mode in 
the spectrum of excitations — gapless means that both 
real and imaginary part of the excitation energy go to 
zero at zero momentum. However, in both cases, the ef- 
fects of pump and decay are such that the real part of 
the Goldstone mode energy is zero also in a finite inter- 
val at small momenta, i.e. the spectrum becomes diffu- 
give^JLL^. This means that a strict application of the Lan- 
dau criterion would lead to a zero critical velocity, where 
quasi-particles can be excited at any value of the fluid 
speed. Therefore, if one would define superfluid proper- 
ties through a strict application of the Landau criterion, 
neither non-resonantly pumped polaritons nor paramet- 
rically driven polaritons in the OPO regime behave as 
superfluids. Nevertheless, in both cases there have been 
evidences for superfluid behaviour. For non-resonantly 
pumped polaritons, it has been recently showrP^ that, 
even though strictly speaking there cannot be superfluid 
behaviour, there are regimes close to equilibrium, where 
the drag force exerted on a small moving defect shows a 
sharp threshold at velocities close to the speed of sound. 
Otherwise, for shorter polariton lifetimes, the threshold- 
like behaviour of the drag force is replaced by a smooth 
crossover. Moreover, metastability of supercurrents in 
non-resonantly pumped microcavities has been theoreti- 
cally demonstrated^. For polaritons in the OPO regime, 
superfluidity has instead been tested through frictionless 
flow of polariton bullet J^, through metastability of quan- 
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turn vortices and persistence of current^- 16 * 17 !. 

The case of coherently driven polaritons in the pump- 
only configuration strongly differs from the two schemes 
previously described: Here, there is no phase freedom 
any longer, because the polariton phase locks to the 
one of the driving pump. As a consequence, the quasi- 
particle excitation spectrum is always gapped — i.e., 
even when the real part of the spectrum energy goes to 
zero at the pump momentum, the corresponding imagi- 
nary part is non zero. For coherently driven polaritons 
in the pump-only configuration, by making use of a lin- 
earised Bogoliubov-like theory, it has been predicted that 
the Landau criterion can instead be satisfied at either 
strong enough pump powers or small flow velocitie s ' 18 1 19 1 , 
regimes where polaritons display a dramatic reduction in 
the intensity of resonant Rayleigh scattering. For values 
of the parameters where instead the spectrum allows the 
excitation of quasi-particles, Cherenkov-like waves were 
predicted^, and recently observed in the density pro- 
file^. Consequently, experiments in this configuration 
have been analysed in terms of the same theoretical de- 
scription which is valid for equilibrium superfluids. 

In this work, we show that, despite the fact coher- 
ently driven polaritons in the pump-only configuration 
do satisfy the Landau criterion at large enough densities 
or small flow velocities, because of the polariton lifetime 
being finite, the fluid always experiences a residual drag 
force even in the limit of asymptotically large densities. 
We show that, only in the limit of perfect microcavities 
(i.e. infinitely long polariton lifetimes), the residual drag 
force at large enough densities goes to zero, recovering 
therefore the equilibrium limit. Otherwise, for finite po- 
lariton lifetimes close to the current experimental values, 



we find that, similarly to the case of incoherently driven 
polaritons, both the drag force exerted on the polari- 
ton fluid by a defect, as well as the height of Chercnkov 
radiation, and the percentage of particles scattered by 
the defect show a smooth crossover rather than a sharp 
threshold-like behaviour which is typical of superfluids 
which obey the Landau criterion. 

The paper is organised as follows: in Sec. [D] we intro- 
duce the model describing polaritons coherently driven 
in the pump-only configuration in presence of a defect 
potential. In Sec. |III| we describe the methods we use for 
our analysis, in particular the numerical algorithm used 
to evaluate quantities such as the drag force, the height 
of Cherenkov radiation, and the percentage of particles 
scattered by the defect which characterise the crossover 
from a superfluid-like to a supersonic-like behaviour. In 
addition, in Sec. |III A} we short ly intr oduce the linearised 
Bogoliubov-like theory of R e f s )iSliiD which will be used 
later in Sec. |IV| in order to compare the results obtained 
with the non-perturbative method with the results ob- 
tained in the linear response approximation. Results are 
discussed in Sec. IV while conclusions (together with a 
discussion of the experimental relevance of our findings) 
are drawn in Sec. |Vl 

II. MODEL 



We describe the dynamics of the resonantly-driven po- 
lariton systeirPH via a Gross-Pitaevskii equation for cou- 
pled cavity and exciton fields ipc,x( r >t)i generalised to 
include the effects of the resonant pumping and decay 
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The single-particle polariton Hamiltonian Hq can be di- 
agonalised in momentum space, 
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by rotating into the lower (LP) and upper polariton (UP) 
basis, 
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where tan2#k = flu/[iOx(k-) — a;c(k)], giving the bare 
lower and upper polariton dispersions: 



(3) 



At zero energy detuning between excitons and photons, 
u;c(k) = wx(k), lower and upper polaritons are an 
equally balanced mixture of exciton and photon, i.e. 
sin 2 0k = cos 2 6>k = 1/2. In contrast, at large values 
of the energy detuning between photons and excitons, 
wc(k) — ^x(k) ^> Qr, i.e. 6^ — > 0, the lower and upper 
polariton respectively coincide with the exciton and the 
photon. 

In the following, we will neglect the exciton dispersion, 
^^(k) = u;x(0), and assume a quadratic dispersion for 
photons, wc(k) = wc(0) + j^j, where the photon mass 
is me = 2 x 10~ 5 mo and uiq is the bare electron mass. 
Through the paper, we will consider the case of zero de- 
tuning at normal incidence, lux(0) = ^c(O), and fix the 
Rabi frequency to fl R = 4.4 meV. The parameters Kx 
and kc are respectively the excitonic and photonic de- 
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cay rates. We will fix these parameters in order to give 
a polariton lifetime, r^p = H/klp, 



«ip(k) = k x cos 2 6> k + k c sin 2 k 



(4) 



close to the experimental values. In addition, we will 
consider the limit of perfect cavities klp —¥ in order to 
recover the equilibrium limit. 

In Eq. ([I]), the exciton-exciton interaction strength gx 
can be set to one by rescaling both fields ipx,C and pump 
strength F p by yOfi/ (2gx)- The cavity field is driven 
by a continuous-wave pump, 



F p (r,t) = Ff,a{r)e 



i(\c p -r—uipt) 



(5) 



with a smoothen top-hat spatial profile with intensity 
/ and full width at half maximum (FWHM) a — 
130 [Am. The pumping laser frequency tu p has been chosen 
0.44 meV blue-detuned above the bare lower polariton 
dispersion at the pump momentum W£p(k p ). The po- 
lariton flow current is determined by the pump momen- 
tum k p , which can be experimentally tuned by changing 
the pumping laser angle of incidence with respect to the 
growth direction ip^: 

ck p = cj L p(k p )sin(</? kp ) . 

Finally, in Eq. ([I]) the potential Vd(r) describes the defect 
acting on the photonic field, over which the polariton 
fluid scatters. Specifically, we consider: 



V d 6(r - r d ) 



(6) 



with r& = 7 /mi and V d = 110 meV. Defects can be 
present naturally in the sample's mirror^. Alternatively, 
defects can been artificially engineered by cither growing 
mesas in one of the mirrors^ or by an additional laseiP? 



III. METHODS 

We numerically solve Eq. |l]) on a 2D grid (256 x 256) 
in a 150 /imxl50 /im box using a 5 th -order adaptive- 
step Runge-Kutta algorithm, and evaluate both exciton 
and photon wave-functions ipx,c{ r i t) in the steady-state 
regime. 

We characterise the crossover from a superfluid-like 
to a supersonic-like regime by evaluating three differ- 
ent quantities. Firstly, we consider the normalised drag 
force^ exherted on the flowing polaritons by the de- 
fect 



J dr\iP c (rW 



dr|Vc(r)| 2 W d (r) . (7) 



When shining the laser pump in the x-direction, k p = 
(k p ,0), the density profile will be symmetric under the 
transformation y i— > -y, i*x,c(x,y) = ipx,c{x, —y), im- 
plying that for defect potentials symmetric under y n- — y 
only the a;-component of the drag force can be non-zero. 




s f 



FIG. 1: (Color online) Drag force in the ^-direction, Fd x 
exerted on the polariton fluid by the defect |6| as a function 
of the velocity ratio c s /vf (see Eqs. (12 1 and (14l), where 
tilp = J dr\ipLp(v)\ 2 /SI is the average polariton density in 
the area il, and for different values of the pump momentum: 



0.7 /jm" (black circles), k p = 0.8 /»m~ (red squares), 
k p = 0.9 fim -1 (green up triangles), and k p = 1 fim^ 1 (blue 
down triangles). The exciton and photon decay rates have 
been fixed to nx = kc = 0.22 meV, while the pumping laser 
frequency uj p has been chosen (as in all the Figures) 0.44 meV 
blue-detuned above the bare LP dispersion cj_Lp(k p ). The 
empty symbols indicate the value of c s /vf above which in 
the linear approximation the Landau condition of Eq. ( 13 1 



cannot be satisfied and therefore where in principle one ex- 
pects the drag force going to zero. Inset: Residual drag 
force at asymptotically large polariton densities ulp (i.e., 
Cs/vf — > oo) as a function of the exciton and photon decay 
rates for k p — 0.7 fim" 1 (black solid line) and k p = 1 /im -1 
(red dashed line). 



Moreover, for step-like defects such as only the val- 
ues of the field |-0c( r )| 2 a t a distance r = contribute 
to the integral Q: 

F dx = r /,y? ,, a / #cos0|Vc(r d ,0)| 2 . (8) 

Therefore the drag force measures the degree of asymme- 
try of the photon (or alternatively the exciton) density 
profiles going from a certain angle <j) < 7r/2 ahead (with 
respect to the fluid flow direction) of the defect, to an 
angle 7r — cf> > ir/2 behind the defect. The asymmetry is 
caused by the scattering of the fluid passing through the 
defect. Plots of the drag force are drawn in Fig. [T] 

Second, we characterise the superfluid-like behaviour 
by the suppression of density modulations aroun d the 
defect known as Cherenkov waves — see, e.g., Refs.B^Ml 
and references ther ein for atomic Bose-Einstein conden- 
sates and Refs.llSllSfor coherently driven polaritons in the 
pump-only configuration. Cherenkov radiation is gener- 
ated in the supersonic regime, when the fluid is passing 
a defect at a velocity higher than the phase velocity of 
the fluid elementary excitations. As mentioned later, a 
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simple analysis of Cherenkov radiation can be carried on 
by making use of perturbative Bogoliubov-like analysis: 
Here, in agreement with the Landau criterion, in the su- 
personic regime, the kinetic energy of the fluid can be 
dissipated radiating Bogoliubov modes, giving rise to per- 
turbations which propagate radially from the defect, with 
a characteristic unperturbed region inside a Cherenkov 
cone which is related to the singularity of the Bogoliubov 
mode dispersion evaluated in the reference frame of the 
moving fluicP^. In our numerical non-perturbative anal- 
ysis, we determine the value of the highest crest of the 
waves, h max , above the mean value of the fluid, h mean , 
giving h^ cr = h max /h meaXi , and analyse how h^ eI changes 
by tuning the fluid velocities and densities. 

Finally, a third way of describing the crossover from a 
superfluid-like to a supersonic-like regime is by evaluating 
the percentage of particles scattered by the defect with 
respect to the total number of particles in the systenPS, 
quantity which we label as S ou t- Both ft,^ er and S ou t are 
plotted in Fig. [5j 



A. Linear response to a weak external potential 

In order to connect our numerical results obtained by 
solving non-perturbatively the dynamics of the Gross- 
Pitaevskii equation ([!]) with the Landau criterion, and 
the spectrum of quasi-particle excitations, similarly to 
what is done for conservative weakly interacting Bose 
gases^, we follow the same perturbative Bogoliubov-like 
analysis firs t intr oduced for resonantly pumped polari- 
tons in RefsP^E^ and expand both exciton and photon 
fields above their mean-field spatially homogeneous and 



stationary states e 
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where we are assuming the pump to have a homogeneous 
profile Ff :<J (r) = f. While the mean-field equations 
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allow to determine the intensity of exciton and photon 
fields in absence of the external potential, fluctuations 
above mean-field need to be introduced in order to eval- 
uate the stability of such a solutions as well as the linear 
response to a weak perturbing external potential. 

The spectrum of the quasi-particle excitations can 
be found by introducing particle-like ux,c and hole-like 
vx.c excitations in both the exciton and photon fluctua- 
tion fields: 



e e ux,c-,k + e e v p ' vx,C;k 



and by solving the eigenvalue equation: 

[(uj + u p )I — L] (u x -M uc;k v x -m wc ; k) T = 0, (10) 
where L is a 4 x 4 matrix given by: 
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Eq. ( 10 1 admits four complex eigenvalues for each k 
which we indicate as w + u„ = LP ± (k),UP ± (k). Note 
that the spectrum obtained this way is already evaluated 
in the polariton flow moving frame. This becomes evi- 
dent in the limit where the pump momentum k p is small 
enough so that the LP dispersion can be approximated 
as parabolic, and the UP dispersion can be neglected. In 
fact, here, one can show^that the spectrum of excitation 
reduces to 



LP ± (k) 



V/ • (k - k p ) - iklp 
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where 



m LP 



(12) 



is the polariton fluid velocity at the pump momentum, 
£k = (k - k p ) 2 /(2m L p), m LP = m c /sin 2 # kp is the 

LP mass, A = lo p — W£p(k p ) — 9lp\vlp\ 2 is the pump 
detuning renormalised by the interaction, and g^p = 
gx cos 4 ■ For more details, we refer the reader to the 
original calculation of Rcf. 19 . What we would like to 
stress here is how to generalise the Landau criterion to 
the complex spectrum of elementary excitations obtained 



from Eq. (10). In particular, in the linear approximation, 
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FIG. 2: (Color online) Residual drag force (in arbitrary units) at asymptotically large polariton densities tllp for k p = 0.7 /im _1 . 
Left panel gives the dependence on the barrier height (in meV) for kx — kc = 0.22 meV. Right panel shows the dependence 
on the decay rates kx and kc (in meV) for two different barriers: Vd = —22 meV (blue dashed line) and Vd = —2.2 meV (red 
dashed line). 



we expect that the perturbation introduced by the defect 
is able to excite Bogoliubov-like quasi-particle states with 
momentum k, when the condition 

K[LP+(k)] -oj p < (13) 

is satisfied. Note that, in the limit where the approxi- 
mation (11 1 is valid and at resonance, A = 0, this recov- 



ers the Landau criterion in its original formulation for a 
conservative system, i.e. close to k p , 5ft[LP + (k)] — ui p ~ 
(c s ± f/)|k — k p | and the critical fluid velocity coincides 
with the sound velocity given by the usuaP' expression: 

(14) 



9hpnLP 
m LP 



where n^p — iV'i'pl 2 is the mean-field polariton density. 
In order to be able to draw an analogy with the equilib- 
rium limit, we will later present our results in terms of 
the ratio c s /vf. 



IV. RESULTS 

We now turn to the non-perturbative numerical analy- 
sis of the Gross-Pitaevskii equation and the analysis 
of the behaviour of the superfluid properties of the reso- 
nantly driven polariton system in the pump-only regime 
when changing the laser pump strength. In particular, 
by fixing the pump momentum k p and energy uj p , we 
evaluate the dependence of the drag force tff\ , the height 
of the Cherenkov waves h^ cr and the percentage of scat- 
tered particles S ou t , on the the steady-state average den- 
sity of polaritons, n^p = J dr|"0£p(r)| 2 /f2, pumped into 
the cavity — here, fi is the circle pumping area we are 
averaging the density over. Rather than as a function of 
nLP 1 we present the results as a function of the ratio be- 
tween sound and polariton fluid velocities, c a /vt, defined 
in Eqs. (pi) and pi}. 
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FIG. 3: (Color online) Perturbation of the steady state photon 
density |^o(r)| 2 (side view) and the y = cut of the profile, 
\ipc{x,y = 0)| 2 (solid black line), close to the defect, for val- 
ues of the pump power asymptotically large (c B /vf — ► oo), 
where the residual drag force of Fig. [T] does not have any ap- 
preciable variation. The pump is shined on the cavity at a 
momentum k p — 1 /im _1 , and the three panels, (a), (b), and 
(c), correspond to increasing values of the polariton lifetime: 
H/kx = h/Kc = 1 ps (a), 3 ps (b) and 120 ps (c). In each 
panel the intensity of the pumping laser has been adjusted 
in order to give the same polariton density. Note that in 
this limit, the defect induces an asymmetric perturbation of 
| t/^c (r) | 2 only ahead and behind the defect, similarly to the 
left bottom panel of Fig. [4] The asymmetry disappears for 
perfect cavities. 



We plot in Fig. [T] the drag force exerted on the polari- 
ton fluid by the defect ^ in the x-direction, Fd x 
as a function of the ratio c s /vf (14 1, for different val- 
ues of the pump wavevector k„ 



v p . We find that the drag 
force decreases fast as a function of the polariton density, 
when c s /vf < 1, while, for finite values of the polari- 
ton lifetime, reaches a finite asymptotic value, a residual 
drag force, at large densities, when c s /i)f ^> 1. While 
for conservative superfluid systems the drag force has a 
threshold-like behaviour and in particular, for perturba- 
tively weak delta-like defects, is finite only for superfluid 
velocities above the speed of souncP, we now observe a 
smooth crossover as a function of c s /vf. In addition, as 



6 





1 V 





LP 


LP + 





-50 50 

X [|un] 



l v 



k [|im ] 



FIG. 4: (Color online) Steady state photon density profile, 
IV'c( r )| 2 (left panels) and corresponding quasi-particle exci- 
tation spectrum ( 10 1 (right panels) evaluated in the homoge- 



neous case for the same system parameters of the left panels. 
Here, k v — 0.9 /im _1 , h/nc = 7 ps, H/kx = 120 ps, and the 
values of the pump power correspond to a state where scat- 
tering is allowed in the linear approximation at c s /vf — 0.54 
(top panels) and no quasi-particle excitation in allowed in the 
linear approximation at c s /vf = 1.49 (bottom panels). Both 
pump strength values considered here corresponds to photon 
densities on the upper branch of the bistability curve. 



shown in the inset of Fig. [Tj we find that the residual 
drag force vanishes only in the limit of perfect microcav- 
ities, Kx = &c 0) i- e -i f° r infinitely long-lived polari- 
tons. Note that this recovers the equilibrium limit, as the 
condition kx = nc ~* automatically requires a pump 
strength / — > 0, as the laser pump strength needed to 
reach a given polariton-density in the cavity decreases 
with the increasing polariton lifetime. 

The two main results to be drawn from Fig. [T] are the 
existence of a cross-over instead of an abrupt transition 
and the appearance of the residual drag force. Since these 
results have been obtained with a practically infinite bar- 
rier (Vd = 110 meV ), one can wonder to what extend 
they remain valid for actual potential barriers which are 
much lower. We have repeated the analysis for differ- 
ent values of Vd, including negative ones, and always a 
cross-over occurs for finite decay rates. More interest- 
ing is the behavior of the residual drag force: it always 
quenches linearly with the decay rates as shown in the 
right panel of Fig. for k p — 0.7 jim^ 1 in two cases with 
Vd = —22 meV and Vd — —2.2 meV. The dependence of 
the residual drag force on Vd is shown in the left panel for 
k p = 0.7 /im -1 and kx = kc — 0.22 meV. Calculations 
have been done with negative values of Vd for numer- 
ical reasons. A clear conclusion is drawn from Fig. [2] 
the residual drag force is non zero as far as both barrier 
height and decay rates are finite. 

The origin of the residual drag force for asymptoti- 
cally large values of c s /vf at finite polariton lifetimes 
can be understood in terms of the asymmetry of the 
perturbation generated by the defect in the (e.g., pho- 
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FIG. 5: (Color online) Height of the Cherenkov waves /ig er 
(black circles) and percentage of scattered particles Sout (red 
triangles) as denned in Sec. |llj plotted as a function of the 
normalized density c s /vf. The parameters are the same as the 
ones of Fig. |4] k v = 0.9 ^m _1 , h/nc = 7 ps, S/kx = 120 ps. 
The empty symbols corresponds to the values c s /vf = 0.54 
and Cs/vf = 1.59 chosen to plot respectively the top and 
bottom panels of Fig. [4] 



tonic) wave- function, \tpc( r ) 1 2 — see Fig. U While, 
as explained later, we observe the disappearance of the 
Cherenkov waves for large enough polariton densities, 
even at asymptotically large polariton densities we al- 
ways observe two small perturbations in the fluid wave- 
function around the defect, one just in front and one just 
behind the defect in the x-direction, i.e., in the polari- 
ton flow direction. We plot in Fig. [3] cut at y = of 
the wave-function around the defect. It is clear from the 
definition ^ that, in this limit, the drag force measures 
the degree of asymmetry of such two perturbations. As 
the perturbation becomes symmetric around the defect 
in the limit of large polariton lifetimes, the drag force 
vanishes. We find that the shorter the polariton life- 
time the smoother is the crossover observed in the drag 
force from supersonic to superfluid behaviour. Instead, 
for high quality cavities, i.e. long polariton lifetimes, 
the transition from supersonic to subsonic behaviour is 
sharper and the residual value of the drag force is smaller. 

In Fig. [T] for every value of the pump momentum 
k p , we also identify the region of values of c s /vt above 
which in the linear approximation the Landau condition 
of Eq. ( 13 1 cannot be satisfied, and therefore where in 



the linear approximation one would not expect any drag 
(such values are plotted as empty symbols). Two ex- 
amples of the linear excitation spectrum are plotted in 
Fig. [4] together with the corresponding profiles evaluated 
beyond the linear approximation — Note that in order to 
do that the pump strength / in the mean-field equations 
of Sec. |III A| needs to be fixed so that to give the same 
average density of polaritons evaluated by solving numer- 
ically Eq. 0, |V>x,cl 2 = / dr\^ x ,c{Y)\ 2 /fi. We find that, 
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gether with the spectrum of linear excitations obtained 
by solving the eigenvalue problem (10) (right panels). 
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FIG. 6: (Color online) Steady state photon density profile, 
|t/>c(r)| 2 (left panels) and corresponding quasi-particle exci- 
tation spectrum (right panels) for k p — 1 /im _1 , S/kc = 
H/kx = 7 ps, and values of the pump power giving c s /vj = 
0.05 (first row), c s /vf = 0.7 (second row), c s /vf — 1.25 (third 
row) and c s /vf — 1.5 (fourth row). Note in that in the third 
panel, while in the linear approximation no scattering and 
therefore no Cherenkov waves are allowed, the full solution 
of the problem still allows scattering, as Cherenkov waves are 
observed in the profile of the corresponding left panel. 



while eventually the Cherenkov waves disappear at large 
enough densities, the persistence of asymmetric pertur- 
bations around the defect, which we ascribe to finite life- 
time effects, contributes to give a finite drag force. This 
is also apparent in Fig. [5j where we plot the height of 
the Cherenkov waves as a function of c s /vf for the same 
system parameters as Fig. |4j Here, it is clear that the 
the Cherenkov waves height is strongly suppressed for 
c s /vf ~ 1, and go to zero for c s /vf 1. 

In addition, in Fig. [5] we plot the percentage of the 
particles scattered by the defect. Here, like for the drag 
force, we find a residual value of the percentage of scat- 
tered particles at asymptotically high polariton densities. 
The difference with the drag force is here that, even for 
perfect cavities when the drag goes to zero, the percent- 
age of scattered particles keeps retaining a residual value 
(see panel (c) of Fig [3]). 

In Fig. \6\ we plot the photon density profiles |i/'c( r )| 
(left panels) for increasing pump power, obtained by solv- 
ing the time dependent Gross-Pitaevskii equation (JlJ , to- 



Left panels show Cherenkov waves evolving smoothly 
from a 'closed' to an 'open' shape till they disappear 
when the subsonic superfluid regime is reached. In par- 
ticular, the angle formed between the waves and the 
propagation direction increases by increasing the polari- 
ton density. Note also that, while qualitatively the phe- 
nomenology of the Cherenkov waves seems to be well de- 
scribed by the linear approximation theory, the transition 
from the supersonic to the subsonic superfluid regime is 
not. In particular, we find a value range of the pump 
power (like the one shown in the third row of Fig. |6j, 
where no scattering and therefore no Cherenkov waves 
are allowed in the linear approximation description, while 
the full numerical solution of the time dependent Gross- 
Pitaevskii equation still allows scattering of polari- 
tons and the associated Cherenkov waves. 



V. CONCLUSIONS AND DISCUSSION 

In this paper we have proposed three different ways to 
analyse the superfluid properties of coherently driven po- 
laritons in the pump-only configuration in presence of a 
defect potential. We have evaluated the drag force when 
the fluid passes the defect, the height of Cherenkov waves, 
and the percentage of particles scattered by the defect. 
By making a comparison with t he linearised Bogoliubov- 
theory introduced in Refs.LUiUiD we have found that, the 
disappearance of the Cherenkov waves, characterising the 
transition from a supersonic to a subsonic superfluid be- 
haviour, is not well described by the linear theory. In ad- 
dition, we have found that non-equilibrium effects which 
go beyond the linear approximation cause a finite residual 
drag force even at asymptotically large polariton densi- 
ties, where the Landau criterion predicts that no elemen- 
tary excitation can be emitted by the defect. Only in 
the limit of infinitely long polariton lifetimes, the resid- 
ual drag force at large enough densities goes to zero, re- 
covering the equilibrium limit. The drag force exerted 
on the polariton fluid by a defect, as well as the height 
of Cherenkov radiation, and the percentage of particles 
scattered by the defect show a smooth crossover rather 
then a sharp threshold-like behaviour which is typical of 
superfluids obeying the Landau criterion. We have char- 
acterised this crossover as a function of the fluid velocity, 
the polariton density and the polariton lifetime. 

The three observables which we here evaluate theoreti- 
cally, can in principle be measured in current state-of-the 
art experiments on semiconductor microcavities. For ex- 
ample, the defect can be carefully engineered by either 
patterning a metal grating on the microcavity top mir- 
ror or growing mesas in one of the mirror^H. This allows 
having a predetermined shape and size of the defect, suit- 
able for a direct comparison with our theoretical analysis. 
Alternatively, the defect can be switched on and off exter- 
nally by an additional laser^l The second scheme would 



allow a direct comparison between the fluid motion in 
presence and absence of the defect. Both schemes are 
within the current experimental reach, and so with this 
work we intend to motivate further experimental inves- 
tigations, which would lead to a better understanding of 
the novel non-equilibrium supcrfluid phenomena in mi- 
crocavities. 
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